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Abstract 

Starting from the new sources of LL- and LR- scalar (slepton) mixings due 
to the R-parity violation, we discuss the structure of lepton-flavor violation 
focusing on the radiative decay of a muon into an electron. Using an optimal 
parametrization, we give the general formulas for the one-loop contributions, 
from which we discuss all combinations of R-parity violating parameters that 
possibly have a substantial contribution to the \x — > e 7 decay width. An exact 
numerical study is performed to obtain explicit bounds on the parameters 
under the present experimental limit. The most interesting one involves a 
combination of bilinear and trilinear couplings. 
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I. INTRODUCTION 



The standard model (SM) is built upon the fact that neutrinos have no right-handed 
components and are massless. Thus, lepton flavors are conserved. So far, experiments 
have not yet observed any lepton-flavor violation (LFV). Nevertheless, the evidence of small 
neutrino masses (implied from oscillation analyses of solar and atmospheric neutrino data 
[ij) implies a small amount of LFV. Many new physics models beyond the SM predict a 
certain amount of LFV, in relation to neutrino mass generation or otherwise. An important 
criterion for a viable model is that it should give an acceptable neutrino mass spectrum while 
staying within the experimental limits of LFV. Among the latter processes, the radiative 
decay of muon, n — > e 7, provides the most stringent bound . The current limit is 

B(fi ->■ ei) < 1.2 x 1(T U . 

A proposal is underway, which can improve this limit by three orders of magnitude Q. 

In the Lagrangian of the minimal supersymmetric standard model (MSSM), apart from 
the soft supersymmetry (SUSY) breaking part, there are no lepton-flavor violating terms 
that would induce lepton-flavor changing processes such as /x — > ej. It is, however, well 
known that soft SUSY breaking parameters can have all types of flavor-changing terms. 
Limits on such parameters from n — > cj have been extensively studied 0. On the other 
hand, flavor-changing soft SUSY breaking parameters are often expected to be suppressed 
by the source of SUSY breaking of a more fundamental theory. A particularly interesting 
example is the gauge mediation models ||. 

Models with "right-handed neutrinos" provide other interesting alternatives with poten- 
tially large LFV 0. It is no surprise that issues of LFV are often tied up with the generation 
of neutrino masses. Within the MSSM, apart from the soft terms, both LFV and neutrino 
masses are indeed forbidden by an ad hoc discrete symmetry — R parity. Note that the soft 
terms by themselves still conserve total lepton number, and hence do not generate neutrino 
masses. 

In the supersymmetric standard model without R-parity imposed, there is however an- 
other important source of LFV. A major part of this comes simply from the R-parity violating 
(RPV) terms in the superpotential. We would like to emphasize that the latter is indepen- 
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dent of SUSY breaking and mediation mechanisms. In fact, the only theoretical means to 
suppress the flavor- changing effects from RPV couplings would be a generic flavor theory, 
or some ad hoc symmetries such as R-parity itself. Hence, it is important to obtain all the 
available experimental bounds on the RPV couplings. We show that the experimental limit 
on the radiative decay of muon, \i — > e 7, provides one of the most stringent constraints on 
the RPV couplings. An improvement of a few orders of magnitude in the experiment may, 
in fact, discover SUSY as well as the LFV, if the model under consideration does explain 
the nature. 

We first explain a new contribution to the Li?-slepton mass mixings involving the bilinear 
fii and trilinear A^ RPV superpotential parameters. Such an important interplay between 
the bilinear and trilinear couplings was first discussed in Refs. in the context of neutrino 
mass. There is also an analogous contribution to the Li?-squark mixings, with implications 
to 1-loop neutrino masses [[J and fermion electric and magnetic dipole moments fT0|-13 . 
With nonzero fi k and A^,,-, the (off-diagonal) Li?-mixing involving the i,j generations is 
proportional to ix% Xkij (sum over k). With no particular arguments or theories to enforce 
Xkij to be diagonal in the last two indices, the LFV proportional to fi* k \ij turns out to 
be quite natural. The combination \x* k X^j conserves the overall lepton number but violates 
lepton flavor. There is another similar contribution to slepton masses in the LL mixing part 
in the form of \x* \iy So far as we know, the implication of such terms in the slepton mass 
matrix is discussed here for the first time. The structure also illustrates very clearly the 
basics of the LFV structure, which comes into fj, — > e 7, for example, in other ways too, as 
we will see below. 

In this article, we presented a systematic analysis of the /i — > e 7 emphasizing the new 
RPV contributions. We work within the framework of the single- VEV parametrization 



(SVP) |i^|T5|1 , which is an optimal choice of flavor basis that helps to guarantee a consistent 
and unambiguous treatment of all kinds of admissible RPV couplings and to maintain a 
simple structure for RPV effects on tree-level mass matrices for all states including scalars 



and fermions [[HJ. Note that under the SVP, the above-mentioned terms with LFV represent 
the full RPV contribution to the corresponding entries in the (charged) slepton mass matrix. 
Obviously, the n* jij combinations contribute equally to the neutral scalar ( "sneutrino" ) 
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masses. There are other interesting RPV contributions that mix the sleptons with the other 
scalars ("Higgses"). Among the latter are the interesting contributions from the soft 
parameters. They come into jx — > e 7 at a different level. 



The present work differs from previous ones on the topic |l7j-j2lj in a few important ways. 
Ours is based on the most generic theory of SUSY without R-parity, i.e., no assumptions on 
the form of R-parity violation, while most of the previous authors worked in various limited 
RPV scenarios. The only exception is the recent paper of Ref. |[20|| , which basically works 
under the same framework as ours. They include the Lf?-slepton mixing contribution, fol- 
lowing the suggestion of Ref. [10(1, where the closely related subject of fermion electric dipole 
moments is discussed. However, the results presented are not as complete and systematic as 
our treatment here. In particular, there is no discussion of the LFV from the bilinear RPV 
parameters only. In addition, their bounds on /4 Xkij have a cos/3 dependence, which we 
disagree with. Other than that, the above new RPV contributions have not been studied 
before. 

The organization of the paper is as follows. In Sec. II, we give the general formulas in 
the basis of mass eigenstates. The focus is on the £~ — > l~ 7 amplitude from 1-loop diagrams 
without colored intermediate states. Though the small deviations of the external physical 
charged lepton mass eigenstates (if 's) from the electroweak states if 's are neglected, full 
attention is given to such deviations for the particles running in the loop when we analyze the 
diagrams. Likewise, exact mass eigenstates are always used in the numerical calculations 
of the 1-loop amplitudes. The formulas can be easily adapted for the radiative r decay 
with simple substitutions. Guided by the formulas and results from our exact numerical 
calculation, we analyze in detail the important pieces of 1-loop contributions in Sec. III. 
Numerical results are presented and discussed in Sec. IV; after which we conclude the article 
in Sec. V. For the readers who are more interested in the major features of the result, rather 
than the sophisticated details on how each RPV parameter comes into play in the process, 
Sec. Ill can simply be skipped. The section also serves as an analytical confirmation of the 
correctness of the numerical results, for the more serious readers. 
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II. CALCULATION IN MASS BASIS 



Here we calculate the decay rate of \i — > through colorless 1-loop diagrams using 
exact mass eigenstates for particles appearing inside the loop. However, as stated above, 
for the external legs to the loop, we use only electroweak states (If s). Each of the latter 
has a minimal deviation from the corresponding L- or f?-handed component of the physical 
charged lepton states as a result of the corresponding small RPV ^ terms . 



A. Model background 



We work in the framework of the generic supersymmetric standard model [15|J2!J , known 
as SUSY without R-parity. Here one could simply give up the notion of R-parity; however, 
the latter terminology is convenient for comparison with the more popular MSSM, as well 
as other limited studies of R-parity violations. We adopt exactly the same notation as in 



Refs. [[Tl],[nj, to which readers are referred for more details. Let us summarize some of the 
useful notation here, starting from the generic superpotential 

W = e ab \^H a u L b a + ylQIH'M + X ajk L a a Qp c k + ~ X^L^l 

(with a,/3 going from to 3, i, j, k from 1 to 3, and e 12 = — e 21 = 1). The essential features 
of the SVP adopted include (pi) = and the identification of L as Hd', related to that is 
the fact that X a = y e . is the diagonal "leptonic Yukawa" couplings, while X ij = for i ^ j. 

(i) We have five (color-singlet) charged fermions from R- and L-handed mass 
eigenstates written, respectively, as ( x+n ) = V T [— iW + , h+, t[, /+, ] T and (x- n ) = 
C7' [— iW~, l~ , /" , Z~ , ] T , where the notation for the electroweak states is quite obvious. So, 
V^Aic U = diag{M-} = diag{M cl , M c2 , m e , m„, m T }. The first two mass eigenvalues are 
the chargino masses. Note that under the SVP, the only RPV parameters going into the 
tree level mass matrix Ai c are the three /ij's. 

(ii) The neutral fermion mass matrix 
Aij^r is written in the basis (— iB, — iW, , W d , Z° , l° 2 , 1° ), where W d = /°. The symmetric, 
but generically non-hermitian, mass matrix is diagonalized by a unitary matrix X such that 
^M-^-X = diagjil^o}, with n — 1 to 4 being the heavy states (neutralinos) and n = 5 to 
7 the physical neutrino states at tree-level. 



(iii) There are 1 + 4 + 3 charged scalars that contribute one unphysical Goldstone state af- 
ter diagonalization. We use the basis { h$, l~, l~, l~, l~, l\\ l%\ It} } to write the mass-squared 
matrix as follows : 



where 



/ Ml MH M* X 
M 2 LH Ml Ml 

•M-RH Mi RL Ml RR 



(1) 



Ml = ml + + Ml cos2/3 - - sin 2 6W 

+ Ml sin 2 /? [1 - sin 2 ^] , 
Ml = m 2 L + m\m L + (/i>^) + Mj cos2/3 

( Ml cos 2 /9 [1 - sin 2 ^] lx3 ^ 



+ sin 2 ^ 



Ux4 



+ 



V 



■Ml = ml + m E ml + Mj cos2/3 [- sin 2 6W 



L 3x3 



(2) 



and 



Ml 



M 2 RH 



(Mir 



\ Ml sin2/3 [1 - sin 2 ^] 



\ 



V2 



(no sum over A;) 







A B / v2 V2 



[A e - /x* tan/3] 



\ V2M W cos/3 



+ 



#2 



/ ^ 



—y* k m k tan/3 



(3) 



Introducing the diagonalizing matrix T> 1 , we have D^Ml'D 1 = diag{ M? m ,m = lto8}. We 
label the unphysical Goldstone mode by m = 1. With relatively small RPV parameters the 
Ml is predominantly diagonal, apart from the mixing between the Higgs bosons (i.e., h u 
and hd = l ) to give the m = 1 mode. The matrix V 1 may then be naturally organized with 
all diagonal entries being of order one and only the 12 and 21 entries being possibly large 
(of order one) among the off-diagonal ones. 
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(iv) The neutral scalar mass terms, in terms of the (1 + 4) complex scalar fields, 



can be written in two parts — a simple (M 



^ )mn ( T m ' 



f> n part, and a Majorana-like part in 



the form | (7W^) mn m (^ n + h.c. As the neutral scalars are originated from chiral doublet 
superfields, the existence of the Majorana-like part is a direct consequence of the electroweak 
symmetry breaking VEVs, hence restricted to the scalars playing the Higgs role only. They 
come from the quartic terms of the Higgs fields in the scalar potential. Using the basis 
( Wj , !° , Z° , l° 2 ,l° 3 ), we have, explicitly, 

sin 2 /? — cos/5 sin/5 lx3 ^ 

1 



Ml 



2 z 



sin 2 /? — cos/3 sin/? 
— cos/? sin/? cos 2 /? 



V 



0, 



0, 



1X3 

o lx3 

0,., 



(4) 



and 



( ^ 



Ml 



V 



/4X + M z 2 cos2/? 



-(B n 



+ M1 



(5) 



m? + (jjCUp) + M| cos2/? 

Note that M.^ here is real (see Appendix A of Ref. ||), while -M 2 ^ does have complex 
entries. Writing the five n 's in terms of their scalar and pseudoscalar parts, the full 10 x 10 
(real and symmetric) mass-squared matrix for the real scalars is then given by 

^ ''VI VI ... 1 

(6) 



Mi 



M 2 SS 



SP 



2 \T 
SP) 



Mi 



where the scalar, pseudoscalar, and mixing parts are 



M 2 SS = MM 2 a ) + M, 



4«t> 



Mp P = Re(A^ 2 1 ) - Ml 



M 2 SP = -MMlt] 



(7) 



respectively. If Im(A^ 2 0t ) vanishes, the scalars and pseudoscalars decouple from one another 
and the unphysical Goldstone mode would be found among the latter. Note that the B a 
entries may also be considered as a kind of LR mixings. 

As a real scalar mass matrix, M 2 could be diagonalized by an orthogonal matrix V s . 
However, we will write V s as if it is just a unitary matrix. This would be useful for illustrating 
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some theoretical features in the discussions below. In fact, it helps sometimes to think about 
the neutral scalars as complex scalars instead of in terms of the scalar and pseudoscalar 
constituents. This is especially true for the "sneutrino" parts. Hence, we write T> sj; Ai? T> s = 
diag{ M? ,rn = ltolO}. Again, it is useful to consider the form of T> s very close to the 
identity matrix, i.e., with all diagonal entries being of order one. The unphysical Goldstone 
mode has, of course, to be found mainly among the first two pseudoscalars. The mode is 
naturally labeled as the m = 6 mass eigenstate here. Similar to the previous case, all the 
off-diagonal entries, except those related to mixing of the Higgs bosons (i.e., the 12, 21, 67, 
and 76 entries) are expected to be relatively small. 

(v) In the diagonalization of the neutral- and charged-scalar mass-squared matrices given 
above, one has to impose all the proper tadpole conditions to get the correct unphysical 
Goldstone modes explicitly. Under the SVP, the tadpole conditions are (see Appendix A of 
Ref. for more details) 

1 



B cot/3 = [m^ + mX + g(<?? + £W ~ %)\ , (8) 
B tan/5 = m 2 o + l^ 2 + ±(g 2 + g 2 )(v 2 - v 2 u ) , (9) 
Bi tan/3 = ir£.+^ (i = 1,2,3). (10) 



The last equation represents the redundancy in parameters explicitly identified under the 
optimal parametrization (SVP) used. Hence, within the formulation, the RPV parameters 
BiS and not exactly independent quantities. Here, the play an interesting 

role very different from the R-parity conserving elements Jtj^.'s in the same matrix. What 
values the fh^.'s might naturally take is quite a technical issue. While the off-diagonal 
mf.^.'s characterize LFV from soft SUSY breaking, the fri^s could be so interpreted only 
when the SUSY breaking mechanism naturally gives vanishing VEV's for the Lj's. Such a 
scenario is indeed possible |23| . In a more generic situation, one can still think about taking 



the low-energy Lagrangian and perform the necessary rotation to recast the model into the 
SVP framework. Such a rotation may generate nonzero f^/s to guarantee the satisfaction 
of Eq . (|10"D in the new basis. Hence, when we switch off the w,f y 's to single out the RPV 
contributions to LFV below, the m 2 's will not be handled in the same fashion. We will 
discuss more the implications of Eq. flTU| ) in relation to fi — > e 7 below. We emphasize here 
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again that the above tadpole conditions are very important and should not be overlooked 
in any discussion of the phenomenology particularly those related to the B^s. 



B. The £j — >• 7 amplitude from (colorless) 1-loop diagrams 



Once we have the above matrices we can express the effective interactions involving an 
external charged lepton with internal particles in terms of the mass eigenstates and the 
elements of diagonalizing matrices. The effective interaction for an external charged lepton, 
taken as an l~ or If state for the L- or f?-handed component here, with exact physical 
charged scalars and neutral fermions inside the loop is given by 



C = g 2 



N L 1 -^+Af R 1 + 75 



2 mm 2 

where |(1 =F 7s) are the L- and Z?-handed projections and 
M R = 



*( x °)</r +h.c. , (li) 



— j= [tan^X^ + X£ n ] T> l ^ i+i)m - ^ X£ n V\ i+ ^ m ^ X* fc+4)n T> l ^ h+5 y m , (12) 

V ^ i/2 i/2 

K L nm = tan9 w X ln T>[ i+5)m - ^ X4nV\ i+2)m + ^ X( i+4 ) n V l 2m -X"( fe+4 ) n V l (h+2)m , (13) 

i/2 i/2 i/2 

with n runs from 1 to 7 and m from 1 to 8. We called this class of contributions neutralino- 
like, which obviously does include the ones with the physical neutralinos among the most 
important parts. The terms in Af^ R are easy to understand. For example, in Af* m , the 
first two terms denote the interactions with L-handed sleptons and the bino and wino, 
respectively. The third term is the interaction with i?-handed sleptons and the higgsino, 
while the fourth term describes the interaction with i?-handed sleptons and the neutrino 
flavor states. Next, we come to the chargino-like contributions. The corresponding effective 
interaction for the external charged lepton If with exact physical neutral scalars and charged 
fermions inside the loop is given by 



C = g 2 

where 



C L 1 -75 C R 1 + 75 



2 ' trim 2 



Vfa) C + h-c. , (14) 



cf nm = -v ln -L [v° i+2)m + iv( i+7)m \ - 1 v {t+2)n -L K _ iV ? m \ 
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92 



1 

V2 



\h+2)m i D(h+7)m] ' 



c l _ a jr/ 2n [P ( S J+2)m + i V* i+7)m ] - j-U (l+2 )n -7= [Am + * V 7m] 



+ 



A, 



klii 



(15) 



(16) 



Recall that T> s is actually real, though we are using V s * notation as if it is not. This is just 
a convention for tracing the LFV structure of the various contributions in our analytical 
discussions below. Here, in fact, the real difference between the T> s * and T> s terms is given 
explicitly by the different signs between the corresponding scalar and pseudoscalar parts. 
Note that the y e . terms in the above expressions can be written together with the A terms 
using the X af3 k notation and the identification of y ei as X u. This common structure between 
L and the Lj's is very useful in our discussions below. 

In applying the above interactions to the process £j(p) — > ^7(9), we can write the 
amplitude as 

1+75 



T = e e* a Ui(p - q) 



m h % a a/3 q 



Ha 



* 1 -Jl + J% 



Uj(j)) 



(17) 



where e* = e*(q) is the polarization four vector of the outgoing photon. The decay rate is 
then simply given by 



ml (\A\\ 2 +\A 



H|2 
2l 



(18) 



It is straightforward to calculate the contributions from 1-loop diagrams with the effective 
interactions of Eqs. ([□]) and ([14]). The result for A\ (Af = A^Il^r) is given by 



AZ, ^em 1 

2 ~ 8tt sin 2 9 w M~ 



-2 



, (M\\ 
M L M F 2 Hf + M R N R 

inm jnm \ ]\/[ I ' tnrri jnm 



where 



a 1 
8tt sin 2 ^ M~f 



Ira , 

f2 



m, 



( M 2 \ m 

mm jnm J I 71 /T2 I %nm l nm rm 

\ lvl s™ I ,n U \ 1V± S 




Mo I .1/7; 

M L M R fJ 



m 



+ C L C 

mm j 



•R* 



jnm 

m 




1 



F 2 (x) = — — — -(1 -6x + 3x 2 + 2x 3 -6x 2 \nx) 
2V ' 6(1 - a;) 4 V ; 

F ^( x ) = n 1 (l-^ 2 + 2x lux) , 
(1 — x) 6 

■ 3 x — 6 x 2 + x 3 + 6 x In x) , 



(2- 



(1 — X) 

F 6 (x) = —^—r{-3 + Ax-x 2 -2 \nx) , 
(1 — x) 6 
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with summations over all physical fermion and scalar mass eigenstates as represented by the 
n and m indices assumed. Apart from the model being different, our background notation 



here, such as the loop integral functions, follows mostly that of Ref. [p4| , to which readers 
are referred for a comparison. 



III. ANALYSIS OF THE VARIOUS CONTRIBUTIONS 

Comparing C L ' R and J\f L ' R , we see that the two types of loop contributions as given in 
Eq.(|l9D do have very similar structures. The first type, corresponding to diagrams with 
a charged fermion and neutral scalar in the loop, are, however, typically larger than their 
£77(2) counterparts, i.e., from diagrams with neutral fermions and charged scalars. Hence, 
we focus our discussions on the C L ' R part, the chargino-like loop. Before going into the 
analysis, it is instructive to introduce the lepton-flavor numbers L e , L^, and L T to the 
superfields as one does to their corresponding components in the SM. Some of the RPV 
parameters would then bear violations of the lepton-flavor numbers. It is obvious that in 
order to have a contribution to [T — > e~ 7, a term must reduce L M and increase L e by exactly 
one unit while leaving L T unchanged. This simple but useful rule serves as a countercheck 
of individual contributions discussed below. 

The following discussions are concerned with sophisticated analytical details on the more 
interesting individual pieces of contributions. The aim is to understand the exact role 
played by each of the RPV parameters in the process, and the strength of each contribution. 
Reading the section would take a bit of effort, which is not necessary for a general reader 
interested mainly in knowing the basic features of the results. While we will refer to some 
analytical results presented here in the discussion of the numerical results in the next section, 
for comparison and confirmation, a general reader may simply skip this section. 



A. Chirality flip inside the loop 

Take A\, for example, among the three terms in Eq. (|T9|) , the one with C\ nm C^ m cor- 
responds to the diagrams with the chirality flip on the internal fermion line. Hence, the 
diagrams with the two chargino states (i.e., n = 1,2) would potentially dominate over the 
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other diagrams with external chirality- flipping mass insertion (m M or m e ). 

We first look at the contributions with a (/i*A) structure. Taking only one A-coupling 
vertex, one can then take the gauge coupling term (first term) from C^ m . We then have the 
real scalar part of the contribution proportional to 



EEKn M- U {k+2)n F 6 ( ^ ) V- V{ h+2)m . ( 20 ) 



n=l m=l 



Here we have dropped the pseudoscalar part just for simplicity in this discussion. This is 
exactly valid when there is no mixings between scalars and pseudoscalars. However, we 
find it more transparent in illustrating the basic features. If the loop function F 6 could 
be factored out from the double summation, we would have a V? M -U(k+2)n summation 
over fermions and a 2^ +2 ) m summation over (real) scalars. However, the fermionic 
sum gives exactly the l^-W + mass term, which is zero. This is well illustrated by the loop 
diagram with electroweak state notation, as given in Fig. 1. Hence, one expects a GIM- 
like cancellation mechanism, here violated only to the extent that F 6 is not universal. The 
violation is, however, quite substantial, as illustrated in the similar case of quark electric 



dipole moments |LT] and by our exact numerical calculation here. The case of degenerate 
charginos is an obvious exception, which, however, needs too large a complex phase for /i to 
be phenomenologically viable. A similar situation goes with the scalar sum — T>f h+2 s = 
8h2 by unitarity. Taking h = 2 in the above expression (|20|) , we have the two dominating 
chargino contributions, the n = 1 and 2 parts, given approximately by 

A, 

& w) ; (21) 

where R R is a 2 x 2 matrix with order one matrix elements and is the i?-handed transformation 



needed to diagonalize the 2x2 MSSM chargino block (see Appendix A of Ref. JTTJ for details.) 
The expected combination ^i* k Afc 21 comes up, with k = 1 and 3 admissible. There is only 
a partial cancellation between the two parts in general. The same situation goes for the 
^inm^2nm P ar t of A 2 , with the combination ^ A£ 12 ( k — 2 and 3 admissible here ) instead. 
An interesting point to note is that the above expression shows no obvious dependence on 
tan/3, a result confirmed by our exact numerical calculation. This is an important issue that 
we will get back to below. 
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Next, we consider the above expression fl2"0|) with h ^ 2. Let us try to estimate the 
individual terms in the T>f h+2 s m summation. From the structure of the neutral scalar 
mass matrix, we see that Ai^s * n particular, is very likely to be predominantly diagonal. 
Off-diagonal elements in the D s matrix then measure the small mixings. The largest terms 
among the different m's are then given bym = 4orm = /i + 2 with a suppression factor of 
magnitude 2% +2 )4 = ^i(h+2) ~ m i ^-m" 1 — • ^his resu ^ is an example of the L-slepton flavor 

m ' L hh mL 22 

mixing from the coupling of the form (X* \Xj mentioned in the Introduction — something that 
is also useful for our discussion below. The flavor mixing is explicitly given in Eq.(|5|), in the 
discussion of the scalar masses, where we explain the notation V s * despite its being real. As 
the RPV mixings are the same in the scalar and pseudoscalar parts and A4p P ), one 

could think in terms of complex scalars here, in which case the complex V s notation would 
be exactly valid. With a first-order difference between mj and mj , the nonuniversal F fi 

L hh L 22 ' u 

function still induce large violation in the unitarity cancellation, essentially between the 
m = 4 and m = h + 2 parts. For the case at hand, the contribution goes with four, instead 
of two, RPV parameters (the two combinations /i 2 \i\ and \i\ Xkhi)', hence, it is very likely to 
be sub dominating. 

There are analogous neutralino-like contributions. However, in the latter case, there is 
also an extra pure gaugino loop, from the first terms in Eqs. (|i~2f ) and (13) with no par- 
allel chargino-like counterparts. This contribution involves V^V^ or V^V^, hence a 
Li?-slepton mixing, with RPV contribution given by fi* k Afc 21 or \i k A£ 12 (for J¥ 2 or respec- 
tively). This is the most intuitive contribution involving the combination of parameters that 
we mentioned in the introduction. Numerically, this contribution is, typically, only similar 
to other neutralino-like terms with one gauge coupling vertex, and hence smaller than the 
dominating chargino-like term discussed above. The reason here is that the larger gauge 
coupling effect is offset by the suppression factor coming from the LZ?-mixing. However, if 
one pushes for a large \/j, \ value while keeping the bino mass M 1 small, the other contri- 
butions, including the whole chargino-like loop contribution could become suppressed, thus 
leaving the pure gauge loop to be the dominating one. 

Replacing the A coupling in the above contributions [cf. expression (^)] with a "leptonic 
Yukawa coupling" (recall y ei = A oii ), we have a 
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5 5 / Ml\ v 

EEK„ M- U 2n F 6 [ 2C V° m V f (22) 

n=lm=l \ iw s m y i/2 

part. Here, the fermionic sum has obviously a dominating chargino contribution without the 
need for RPV parameters. However, we have no choice but the fixed scalar mixings given 
by V^V^. As discussed above, the latter involves /i* fj, 2 - With the small y ei coupling, it 
gives a strong suppression. Of course the "electron Yukawa" y ei above can be exchanged for 
the larger "muon Yukawa" y e2 in the corresponding piece in A% instead. The latter give the 
same X^X> 3 s m scalar mixing factor. This contribution is depicted in Fig. 2. 

From Eqs. (|T5|) and (Jl6|), there is one more similar but independent term. This has 
expression ( |2"2"D modified with a simple switching of the explicit 2 and 3 indices, and a sign 
flip. Here we write down explicitly the corresponding term from A 2 instead : 

5 5 /Ml \ y 

E E -Vm M x - n ir 4n fJ v° m vz y f . (23) 

n=lm=l V Sm / g ' 2 

From the above discussion, it is easy to see that the dominating part with the charginos 
(n = 1 and 2) gives a /i 2 dependence through U^ n . Naively, we expect a /i* from the scalar 
mixing part. The mixing of the corresponding complex scalars, however, involves a 32- 
entry of the -M^a matrix. From the result given in the Sec. IIA, we see that the mixing 
does have a n*. But it comes in the combination rh? + //* // The latter is, by the tadpole 
equation (|10D, B\ tan/?. Detailed analysis of the full scalar mass matrix, with the first-order 
mixings among the Higgs states, does not change the tan/3 dependence of the result. 
Hence, the LFV structure here is given by the RPV parameter combination B\ /x 2 , which is 
related to the /x* /j, 2 combination through the tadpole equation. Given in terms of B^, the 
contribution has an explicit tan/3 dependence. In general, one expects this contribution to 
be of comparable magnitude to the previous one. 

The last contribution we want to discuss here is the one with two A couplings. Extracting 
the part similar to expression (EO), we have 



/ 5 I Mi \ x h , 2k , A hl 

E E V Tfc'+2)n M x -U{k+2)n FA — ^ [V'y +2)m + z£> ( S fo , +7) J [V° h+2)m + ll^+^j — f - 1 

m n =l \ 1V1 S m / 92 92 



(24) 

Note that we have written the expression very differently here. Both the scalar and pseu- 
doscalar parts are explicitly shown. The sum over m, scalar mass eigenstates, goes over all 
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nine physical states (the J2' m notation means the unphysical Goldstone mode is omitted). 
The interesting point here is that this contribution involves a nontrivial interplay between 
two parts. A careful reading of the X* h , 2k , term in Eq.(^) would appreciate that the scalar 
mixing matrix comes in as T> s * instead, or, equivalently, the scalar and pseudoscalar parts 
come in with the "wrong" sign. At the limit of degenerate but unmixed scalars and pseu- 
doscalars, the "right" signs (or with V s , as in the case of the first gauge coupling term) would 
give identical contributions from the two parts; while the "wrong" signs would give an exact 
cancellation instead. To put it in another way, a Majorana-like scalar mass insertion along 
the scalar line is needed to have the contribution nonzero, as depicted in Fig. 3. There 
have been some discussions on the Majorana-like scalar mass terms for the "sneutrinos" 
(or 1° states to be exact) as complex scalars, and some of the resulting phenomenological 



implications are studied in recent literature [25,26]. Such mass terms appear under the SVP 



through LFV from the soft Bi parameters. An illustration of the Majorana-like scalar mass 
terms is shown explicitly in Fig. 3b. It can be easily seen that the contribution under discus- 
sion here actually involves a minimum of four RPV parameters, By B\ and XyiV Xkhi, and 
it is achieved by taking k = k'. An interesting point, at least from the theoretical point of 
view, is that the contribution could have a different flavor structure from the more familiar 



two-A loop diagrams with the chirality flip on the external muon line [ 19 1 , and a m T instead 
of rrifj, dependence too. For instance, we could have a new contribution from the A123 A 311 
combination. 

Finally, we want to remark that the remaining contributions with internal chirality flip 
can be analyzed again by using the extended flavor structure with the "leptonic Yukawa 
coupling" y e . identified as X oii , as discussed above. In fact, Fig. 3a is given with the generic 
A Q(3 fc notation and hence applicable to this latter case. One can obtain contributions de- 
pending only on a B^ and A^i combination, for example. Explicitly, we have expression 
(El) modified to 



' 5 ( M2 -\ v -\ h 

E E Vtn K-U, n Felj^j [V° m + W°J [V { l +2)rn + iV( h+7)m ] ^ — ^ , (25) 

with which the fermionic sum suggests a major contribution from n — 4, i.e., the muon itself 
with the m M dependence, and the scalar sum gives the dominating contribution proportional 
to B^. However, this contribution then has two factors of "muon Yukawa" (y e2 ) suppression. 
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Hence, it is not a very important contribution. Note that the Majorana-like scalar mass 
required can be interpreted as in the same form as given by Fig. 3b, with a Bq and a B^, 
except that the Bq insertion, in general, has no suppression; the Bq is actually not necessary 
to complete the corresponding diagram in Fig. 3a in this case though, due to the existence 
of the I/ (or, equivalently 1° = h d ) VEV. 



B. Chirality flip on the external muon line 

Contributions from the terms with the chirality flip on the external muon line is expected 
to be suppressed by a y e2 (~ 10~ 3 ) factor. However, we have seen that apart from the /x* A- 
type combination, other combinations of RPV parameters that come into the diagrams with 
an internal chirality flip do have extra suppressions. Hence, it is still of interest to see if there 
are important terms with an external chirality flip that are comparable to those discussed 
in the previous section. Therefore, we mainly look for terms depending on only two RPV 
parameters without further extra suppressions, apart from the one dictated by the chirality 
flip. 

For the chargino-like diagrams, for instance, our candidates here are from the C^ nm C^* m 
part of A\ and the C£ m C^* m part of A\. The latter looks more interesting. First of all, there 
is a pure gauge loop, as shown in Fig. 4a. Compared with the above, we have the analogous 
expression 

5 5 ( M 2 \ 

™,EE - K, Vm F 5 [ D| m . (26) 

n=lm=l \ 1VI S m ) 

The scalar mixing part gives the /x* /i 2 LFV structure. This looks comparable with expression 



The only other term from C^ nm C^* m without further y e . suppression is the usual 2-A term 
discussed much earlier in the literature, together with similar (colored) 2-A' diagrams JTU . 
We have an expression given by 

5 5 / Ml \ x h/2k , x* 

m ^EE -V ffc'+2)n V (k+2)n F 5 \ — ^ ^ +2 )m V (h+2)m ~ f ^ . ( 27 ) 

n=im=l \ m s m ) 9i 

requiring further h = h' and k = k'. Note that all off-diagonal matrix elements of the 
form V(fc+2)n are very small, those RPV ones (n — 1 or 2) in particular further contain a 
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y e . suppression (see Appendix A of Ref. [fLl]]). Finally, with the Cf nm C£* m part of A\, we are 
obviously left only with the other 2-A terms 

The neutralino-like loop contributions are mostly straight forward analogs, except for the 
contributions with only one gauge coupling. Unlike the chargino-like case discussed above, 
there may be no extra y e . suppression, from a V element or otherwise. However, an extra 
Li?-mixing is involved, giving another suppression, as illustrated in Fig. 4b. 

The only relatively important type of contributions involving the combination such as 
A/i2i has two factors of y e2 suppression. If one is interesting in the dominating contribution 
of this type, an external chirality flip term with only one extra y e2 suppression may have to 
be considered. In fact, such a contribution does exist. The C, L Ct* part of Ah has a piece 

' lnm 2nm 1 Z J- 

of the form 

5 5 (Ml\ x _ u 

™,EE '"In U, n fJ V( h+2)m VZ — - -p , (28) 

n=\m=\ \ IV1 S m J Ml 

where we have h = 3 or 1. This is expected to be of similar strength to the one in expression 
(p5|) discussed above. Recall that the scalar mixings give a tan/3 dependence. 

IV. EXACT NUMERICAL CALCULATIONS AND RESULTS 

In this section, we present the results we obtained by a careful numerical implementation 
of our /i — > e7 formulas, with explicit numerical diagonalization of all the mass matrices 
involved. We isolate various major contributions by singling out each of the corresponding 
RPV parameter combinations as the only nonvanishing ones. The soft SUSY breaking 
contributions to R-parity conserving slepton mixings are set to zero ( i.e., mf, rhj, and A E 
are set to be diagonal ). Though we have used only real numbers for all input parameters 
in the numerical results presented, our discussions apply to the general case of complex 
parameters, as does the analysis given above. A basic set of typical values chosen for the 
input parameters are given in Table |. We used this set of inputs unless otherwise specified 
in the results below. At the end, we also show the effects of varying these input parameters. 

At the beginning of the previous section, we introduced a useful rule in terms of lepton- 
flavor number violation counts; i.e., an admissible contribution has to include RPV pa- 
rameters such that it reduces and increases L e by exactly one unit while leaving L T 
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untouched. Moreover, 11 — > e 7 is an R-parity even process. The rule alone can be used to 
identify a few interesting cases of two RPV parameter combinations. We first concentrate 
on the superpotential parameters. Using the bilinear parameters, we have only the fx* /i 2 
combination. With the trilinear parameters, there are various 2-A and 2-A' combinations. 
Such contributions to ll — > e-j are well studied [B5|. It is easy to see that as the A' couplings 
involve quarks and squarks while the others do not, there is no combination of a single A' 
with a coupling of the other types contributing to 1-loop \i — ► e 7. We skip the quark-squark 
loop contributions here. The 2-A loops, however, is an integral part of our (colorless loop) 
formulas. We will discuss such contributions also for completeness and for easy comparison 
with previous works. Most interestingly, however, are the combinations involving RPV pa- 
rameters of both the bilinear and trilinear types. They are ii\ Afc 2i with k being 3 or 1, and 
/ifc \* kl2 with k being 3 or 2. 

Next, we add into consideration the RPV parameters from soft SUSY breaking. The 
trilinear A x terms have no role to play here, because under the SVP such terms give only 
three-scalar interaction vertices, which obviously cannot be incorporated into any 1-loop 
fi — > e7 diagram. This leaves the B^s and rn^.'s. However, as commented above in the last 
part of Sec. IIA, we have the important tadpole conditions in Eq. ([TDD relating these two sets 
of parameters to the fx^s. In the literature, fiiS and B^s are usually discussed as independent 
parameters while the WjL's and the tadpole conditions are often overlooked. Recall that 
imposing the correct tadpole equations is crucial in getting the correct physical scalars. If 
one is in favor of setting w.^.'s to zero, a Bi value is then fixed by the corresponding /ij, 
and vice versa. However, in order to highlight some of the analytical features discussed and 
to enable an easy comparison with previous studies of other phenomenological implications 
of the Bi in the literature, we will first discuss the //j's and £?j's separately as if they are 
independent. We will first single out a contribution involving a Lii by setting the Bi to 
zero and tuning rh^. to satisfy the tadpole condition before solving for the spectra of the 
physical scalars. Likewise, we will tune a /Zj to zero to single out a Bi effect. An interested 
reader can put the two otherwise related contributions, for example a //* A 321 and a B^ A 321 
contribution, together by imposing the tadpole equation and picking whatever m^. value 
deemed appropriate. We will also show a case of m? . = for the /1* fi 2 contribution at the 
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end. Please be reminded that in the tadpole equation, Bi goes with a tan/3 factor. For the 
same /Xj, will then be suppressed by tan/?. 

It is easy to see that under the strategy discussed, the additional RPV parameter combi- 
nations of interest are given by B\ /i 2 , /i* B 2 , and B\ B 2 , together with B\ Afc 21 and Bk A£ 12 . 
These complete our list. Now, we go into each of these combinations of RPV parameters. 
Readers interested in more analytical details are urged to compare our discussions below 
with those presented in the previous section, the cross references to which are given inside 
square brackets ([ ]'s). 



A. The /x£ Afc 21 or A£ 12 contributions. 

This is the most interesting case, because it involves both the bilinear and trilinear RPV 
couplings and also it is likely to give a much larger branching ratio. Later, we will see that 
the bounds obtained on such parameter combinations are comparable to what one could 
obtain by imposing a sub-eV bound on all neutrino mass contributions (see, for example, 
Refs [|]J||27|)- Without loss of generality, we take the /i* A 321 combination for illustration. 
The dominant contribution comes from the last term of Eq. (|19"D as shown in Fig. 1 [also 
discussed in expression (EO)]. This is confirmed by our exact numerical calculation. Another 
interesting contribution comes from the pure neutral gaugino with /j* A 321 coming in through 
the Li?-slepton mixing. There are also the neutralino contributions without LZ?-slepton 
mixing, the exact analog of the chargino ones. In a generic region of the parameter space, 
the chargino-like loop result typically dominates over the neutralino-like loop result. We 
plot contours of the resulting branching ratio as a function of (real) /i 3 and A 32 i in Fig. |5|. 
The present experimental limit is also shown and the allowed region at 90% C.L. is shaded. 
The contours for the other three combinations of RPV parameters, each taken alone, are 
essentially the same. Recall that the corresponding dominating contribution for the two 
//fcA£ 12 combinations comes in via the Cf^C^™ P ar ^ °f ^2 instead. The 90% C.L. upper 
limit on |/^Afc 21 | or |/i&A£ 12 | (normalized by |/i | = 100 GeV) is given by 

I/A3 ^321 I 1/^1 -^121 I |/Aj Agiaj \Jh ^212! < 1 5 X 10~ 7 (29) 

I /A) I I /A> I I /A) I I /A) I 



As to be explicitly illustrated below [and discussed with expression (|2lD], this kind of 
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contribution is insensitive to the tan/3. Recall that imposing a sub-eV bound for all neu- 
trino mass terms obtainable from any RPV parameters (see, for example, Ref. [2"7|]) gives 
|M cos/3 < 10~ 6 . A simple estimate of a 1-loop neutrino mass diagram from two A- type 
couplings give the corresponding bounds |A 121 A 212 | , |A 321 A 312 | < 0.015, and |A 212 | < 0.008. 
These bounds are not very strong at all as factors of and m e are involved in the neu- 
trino mass diagrams. Moreover, A 212 is the only parameter that is capable of giving rise 
to a 1-loop diagram just on its own, hence with a bound on itself alone. The other four 
A couplings are otherwise bounded by 0.04 ||28|| , from charged current processes and r de- 
cays. The Hi parameters typically have no strong bound apart from the one due to neutrino 



masses (see Ref. [13 for details), which is suggested but not mandated from the result of 
the super-Kamiokande experiment |29| . Hence, we can see that the bound we obtained here 
from fi — > e7 is very important, especially in the large tan/3 region where the neutrino mass 
bound is weakened. Further improvement in the // — > experiment is capable of giving the 
best bound on the RPV parameters, or discovering the signal of th R-parity violation. 



B. The fi* /x 2 contribution. 

The contribution is dominated by the Cf nm C2^ m term of A 2 as depicted explicitly in 
Fig. 2 [cf. expression fl22"[) 1. The result obviously contains a "muon Yukawa" suppression. 
The same story goes for the contribution from the C^ nm C£* m part of A 2 , depicted in Fig. 4a 
[cf. expression (p6|)] . The latter is typically numerically smaller. We show contours of 
B(n — > ej) in the real (z^,/^) plane in Fig. The 90% C.L. upper bound on h* /x 2 , 
normalized by |/x | 2 (|/x | = 100 GeV here), is 

l/i * / " 2 ' <0.53x lO" 4 . (30) 
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Note that the bound weakens roughly by a factor of ^ compared with , as expected. 

Though the bound looks weak compared with the sub-eV neutrino mass bound discussed 
above, it is still significant when compared with most of the other bounds on the /x x and /j, 2 
parameters, especially in the region of large tan/? [[14] . 
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C. The contributions from two A-type couplings. 



Here we have naively two classes of contributions, a class of A* A diagrams and a class of 
A A or A* A* diagrams. The first class needs a chirality flip on the external muon line. We 
have a A* 3fc X 32 k combination (for k = 1,2,3) from the C£ m C^ m part of A% [cf. expression 
(PT|)1, and a A* J2 A^i combination (for ij = 12,13,23) from the C^ nm C^* m part of A L 2 . The 
first group must have a running in the loop, while the second admits both a /" and a l~ 
and hence gives a twice stronger result (with roughly degenerate L and R sleptons). The 
bounds obtained on appropriate combinations of AA* are given in Table [TT], in which case we 
have deliberately used tan/3 = 1 (instead of 10, in order keep the physical L- and i?-handed 
slepton masses at about 100 GeV) to show the exact agreement with Ref. fl9 |. Note again 
the numerical bounds are roughly a factor of 10~ 3 weaker than the fi* A-type bound, as a 
result of the m M factor. 

On the other hand, a AA or A* A* diagram requires no chirality flip outside the loop. In 
the AA case, for example, we have to pick a \h>2k> to get the required reduction in L M and 
a Xhki to get the increase in L e . It is easy to see then that it is impossible to choose a 
combination of AA such that it does not cause further changes in any of the lepton-flavor 
numbers, unless more RPV parameters are involved. A term of the latter case has been 
discussed analytically [cf. expression Q23j)]. The same situation holds for A* A*. We will not 
further investigate such contributions here. 



D. The contributions involving the Bi parameters 

We have discussed above the implication of the tadpole equations relating the £?j's to 
the /Vs. Our numerical strategy isolates terms explicitly proportional to a Bi or a /^. 
Here we discuss the contributions involving the -Bj's. First, there are combinations B% Afc 2 i 
and Bk A£ 12 . For B% Afai an illustrative contribution has been discussed analytically in the 
previous section [cf. expressions ( p5|) and (p8|)]. The bounds obtained on these combinations 
are shown in Table O. Contours of B(fi — > e 7) in the real (-B3, A 32 i) plane are shown in Fig. |^. 
Understanding the result analytically is more complicated here. Our analysis does suggest 
more than one factor of Yukawa suppression, hence the weakness of the bound obtained. 
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Numerically, the part with external chirality flip from a C^ nm C^* m contribution [cf. expression 
(f2"5|)] can dominate over the part with the chirality flip inside the loop [cf. expression (P5|)] 
depending on tan (3. The overall bound looks stronger than the analytical estimate. This 
is to be explained by the larger loop function F$ with the much lighter muon propagator 
compared with that of a chargino, and the tan/? dependence. However, the relation between 
a B k and a fi k and the weakness of the present result compared to the fi* k \ k2 i or fj> k X* kl2 
result suggests that the B* A- or BX*- type contribution is really of less significance. 

Next, we have the B^ \i 2 combination. Recall that according to our strategy, this is probed 
with fi ± and B 2 set to zero. It clearly has a "muon Yukawa" suppression [cf. expression 
(PH)]. Numerically, the bound is similar to that on /j,* /i 2 . The number shown in Table |Tl 
actually looks better than the corresponding ^\ // 2 number. However, B\ comes in with a 
tan/3 dependence. Therefore, 1^4 should be compared with ^r^pr^ hence explaining the 
difference. Again, we give contours of B(/j, — > ej) in the real (f?i,/x 2 ) plane, in Fig. |8]. 

The analogous contribution of the type coming from /i* B 2 has a "electron Yukawa" 
suppression. Numerically, it is confirmed to be much smaller. Likewise, we find no important 
contribution from an explicit B\ B 2 combination. In fact, one obviously cannot make a 
simple fi — > e 7 diagram with a B\ and a B 2 RPV insertions. 



E. Parameter Variations 



In this section, we illustrate the effects of varying the input parameters on the bounds, 
using Afcai| and /i 2 | as examples. The results are summarized in Table JTJ. Basically, 
the effects of varying the mass parameters (fi , M 1 , rh£ , mj) reflect on what particles are 
inside the loop of the dominant diagrams. In the case of \/a* A 321 |, the dominant diagram 
involves mainly the / 2 , while the l/i* /x 2 | case involves the mixing between / 2 and 1°. Therefore, 
varying m| does not have much effect on the bounds while varying the corresponding entries 
in ml changes the bounds significantly. Obviously, a large (relevant) scalar mass suppresses 
the fi — > e7 amplitude and thus weakens the bounds [see parts (iii) and (iv)]. 

Increasing fi and M 1 = | M 2 essentially increases the chargino and neutralino masses [see 
parts (i) and (ii)]. Here, the variation of the bound is more complicated. For illustration, we 
show more details of the fi variations in part (i). Taking A 321 | as an example, the bound 

22 



is most stringent for small |/i |. However, for negative /i , apart from the general weakening 
trend as |/i | increases, there is an extra structure in the dominating chargino contribution, 
namely, there is a dip at /i = — M, tan/?, where it essentially vanishes. This is a special 
feature of the type of RPV contribution also observed in the similar contribution to neutron 
electric dipole moment [O. For the relatively large value of tan/3, used here, however, this 
is already well inside the large \y, Q \ region where the pure gauge loop contribution becomes 
dominant. The latter case generally happens when the bino mass M 1 is small relative to 
\(j, \. Note that the pure gauge loop contribution is independent of fi . Apart from the dip 
mentioned, the dominant chargino contribution does decrease with increasing \(J, \, as shown 
in column 1 of Table [TIT]. 

As for the \fj,* [a 2 \ case, the variation of the bound with \/j, \ is more complicated. Unlike 
the previous case, the contribution has, analytically, an explicit chargino mass dependence 
[as shown by comparing expression (p2"|) with expression (pT|)1 . Hence, we do not expect a 
simple weakening of the result as |/i | increases. Another important issue here is that there is 
more than one important piece of contribution of this type in most regions of the parameter 
space. For example, if we stick to setting Bi = B 2 = as we do to obtain the numbers 
given in column 2 of Table [TTT| , we still have the pieces corresponding to Fig. 2 and Fig. 4a 
[cf. expressions ( ^2j) and in interference with one another. Besides, at very large |/i | 
(and relatively small M x ) the pure gauge-gauge term of the neutralino-like contribution is 
increasing and dominates over the chargino-like contribution. 



In addition, we have shown in column 3 of the Table |T|, an interesting case with . = 
instead. Then the n\ yu 2 - and Bl /z 2 -type contributions exist simultaneously. As given in 
Eq. (|IUD , we then have Bi tan/? = /i* /u 1; and the contribution of Bl fj, 2 type such as that 
given by expression (|23|) comes along and interferes with that of the /i* /i 2 . Our result in 
Table |H| suggests that the Bl /i 2 type contribution is smaller than the /x* /i 2 type but are of 
the same order; a destructive interference of the two parts then weakens the overall result. 
Note that our analytical discussions in the previous section [cf. expressions ( p2j) and (p3|) 1 
indicate that the terms come in with a different sign. With larger |/x |, the bounds shown in 
column 3 of the table are much weaker, because of a stronger cancellation as the Bl /i 2 term 
increases. This is mainly a result of the increase in the Bi value (for fixed /uj, however, 
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the contribution from the term has more complicated dependence on the \/i \. Similarly, 
the same increase in cancellation with larger , as shown in the table can be understood 
from comparing the two terms. The results shown in column 3, hence, further illustrate the 
importance of the tadpole condition given by Eq.(|i~0|) emphasized throughout the paper. 

Finally, we comment briefly on the tan/? dependence of the results, also illustrated in 
part (v) of Table JT|. From the table we can see that varying tan/3 has only a little effect 
on \fj,* A 3 2i| but a rather significant effect on \fj,* fj, 2 \. The lack of sensitivity to tan/3 in the 



former case has been suggested in our analytical discussion [cf. expression (2"I)] and further 
confirmed over the range of the tan/3 value. In the latter case, the numerical result shows that 
the bound is strengthened by a factor of cos/3. This simply illustrates the ^3 dependence 
of the Yukawa coupling y e2 . 

V. CONCLUSION 

In this paper, we have performed a comprehensive study on the radiative decay of muon 
(/i — > 67) in the framework of the generic supersymmetric standard model (without R 
parity). We have identified a few combinations of a minimal number of RPV couplings 
contributing to the decay. Among them the most interesting are /4 Afc 21 and Hk \% 12 - The 
upper bound on the combinations obtained from the experimental limit on /x — > e 7 is 

<l.5xl0- 7 , 

which is as stringent as the ones that can be obtained from the constraint of sub-eV neutrino 
masses. Note that different combinations of RPV parameters are involved in the generation 
of neutrino masses though. Furthermore, our result, in contrast to a similar result given in 
Ref. [20|, has little sensitivity to tan/3. 



Another combination, /i* /i 2 , contributing to [L — > involves the LL-slepton mixing. 
This contribution is identified for the first time. The upper bound on the combination 
obtained from the experimental limit on fi — > e'j could be important, especially in the 
region of large tan/3. We have also discussed the related role of the soft SUSY breaking Bi 
parameters in the process. 

Before closing we summarize the following important points : 
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(i) The combinations n* k Afc 21 and /i k \* kl2 participate directly in the Lf?-slepton mixings, 
while n* Hj participate in the LL-slepton mixings. The combinations also highlight the major 
/j, — > e7 contributions. Under our formulation (SVP), the only RPV soft SUSY breaking 
parameters that contribute to (tree-level) slepton masses are the Sj's and the rh^s. 

(ii) In relation to the RPV parameters (under SVP), the tadpole equations say that 
{fii, Bi, m^.} for each % are not independent. This fact is often overlooked in the literature. 
In our study, we single out the contribution from fa or Bi with a matching nonzero m^.. 
We have also illustrated a case with m^. = 0, and thus combined the effects of both /ij and 
Bi. This case shows that there could be strong cancellations between them. Hence, the 
overall contribution from the bilinear RPV parameters only is tied up with how the tadpole 
equations among {fa, B iy m^.j-'s are chosen to be satisfied. 
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Table captions : 

Table I — Basic input SUSY parameters for the numerical results presented. These values 
are adopted unless otherwise specified. Note that m^ o corresponds to the soft mass square 
for Hd, in the MSSM language. Moreover, soft masses for H u and Bq (~ MSSM soft SUSY 
breaking B term) are not used as inputs, but are determined from rn^ and fi (~ MSSM \i 
term) through the tadpole equations for correct electroweak symmetry breaking. 

Table II — Summary of bounds on various combinations of two R-parity violating param- 
eters, normalized by |/x | = 100 GeV where appropriate, due to the experimental limit 
B([j, — > e'y) < 1.2 X 1CT 11 at 90% CL. The input parameters are as in Table [I], except for 
the limits on |AA*| combinations in which case we have used tan/3 = 1 (as explained in the 
text). 

Table III — Effects of parameter variations of interest, on the bounds of A 32 i | ■ (100 GeV) -1 
and [i 2 \ • (100 GeV) -2 . Note that the fixed mass scale of 100 GeV is used for normalization 
to extract numerical bounds. The tadpole condition is chosen with B^ = in the second 
column of the bounds while m? = is chosen in the last column. 
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Figure captions : 

Fig. 1 — The R-parity violating chargino-like loop diagram. 
Fig. 2 — Diagram with charged gaugino/higgsino mixing. 

Fig. 3 — A diagram involving Majorana-like scalar mass insertion. 3a/ The diagram. 3b/ 
The seesaw origin of Majorana-like scalar masses for the "sneutrinos" explicitly illustrated. 
Fig. 4 — Illustrative diagrams with chirality flip on the external muon line. 4a/. A charged 
gaugino loop. 4b/ A g — A neutralino-like loop. 

Fig. 5 — Contours of B(/i — > ej) in the (real) plane of (/x 3 ,A 32 i). The 90% C.L. allowed 
region is shaded. 

Fig. 6 — Contours of B(/j, — > e^y) in the (real) plane of (/ii,/i 2 )- The 90% C.L. allowed 
region is shaded. Note that the approximation we used for the external lepton lines is less 
applicable at the right and top ends of the plot where the result should be read with caution. 
Fig. 7 — Contours of B(/j, — > e^y) in the (real) plane of (S 3 ,A 32 i). The 90% C.L. allowed 
region is shaded. 

Fig. 8 — Contours of B(jj, — > ej) in the (real) plane of (Bi,fi 2 ). The 90% C.L. allowed 
region is shaded. 
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TABLES 



TABLE I. Basic input SUSY parameters for the numerical results presented. These values 
are adopted unless otherwise specified. Note that ml 00 corresponds to the soft mass square for Hd, 
in the MSSM language. Moreover, soft masses for H u and Bq (~ MSSM soft SUSY breaking B 
term) are not used as inputs, but are determined from m? Q and /jl (~ MSSM fj, term) through the 
tadpole equations for correct electroweak symmetry breaking. 



M 1 (GeV) 


M 2 (GeV) 


/i (GeV) 


tan/3 


100 


200 


100 


10 


ml (10 4 GeV 2 ) 


ml (10 4 GeV 2 ) 


A e (GeV) 





diag{2, 1,1,1} diag{l,l,l} 100 



TABLE II. Summary of bounds on various combinations of two R-parity violating pa- 
rameters, normalized by \fj, \ = 100 GeV where appropriate, due to the experimental limit 
B(/j, — > cy) < 1.2 x 10 -11 at 90% C.L. The input parameters are as in Table |, except for the limits 
on | A A* | combinations in which case we have used tan/3 = 1 (as explained in the text). 



1/4 


A321I 


1/4 A121I 


1^3 As 12 l 


nr Im2A* 12 | 


< 1.5 x 10~ 7 




Mol ' 


l«)| ' 


ImoI 


' or ImoI 


KH 








< 0.53 x 10~ 4 


1^321 A* 31 | 


IX X* 1 

) l A 322^i 3 2| 


, or 


^323^*33 


< 2.2 x 10" 4 


1 X* X 1 

l A 132 A 13l| 


1 X* X 1 
! l A l22 A 12l| 


, or 


^232^231 


< 1.1 x 10" 4 


|B 3 *A 32 i| 


\B*X 131 \ 
IwP ' 


1^3 


1 nr \B2\* 211 \ 


< 2.0 x 10~ 3 


Iwl 3 








< 1.1 x 10~ 5 
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TABLE III. Effects of parameter variations of interest, on the bounds of |/ig A 321 | • (100 GeV) -1 
and fi 2 \ • (100 GeV) -2 . Note that the fixed mass scale of 100 GeV is used for normalization to 
extract numerical bounds. The tadpole condition is chosen with Bi = in the second column of 
the bounds while m? = is chosen in the last column. 

Parameter changes Normalized numerical bounds 

Ijfa A 32i| [Afi£a| Im*M2| 

(100 GeV) (100 GeV) 2 (100 GeV) 2 











(with Bi = 0) 


(with ml . = 0) 


Original inputs of Table | 


< 1.5 x 10" 


-7 


< 0.53 x 10~ 4 


< 2.3 x 10~ 4 


(i) Mo = 


500 GeV 


< 10 x 10 _ 


7 


< 0.80 x 10~ 4 


< 27 x 10~ 4 


Mo = 


250 GeV 


< 5.8 x 10" 


-7 


< 1.1 x 10~ 4 


< 14 x 10~ 4 


Mo = 


-100 GeV 


< 2.0 x 10" 


-7 


< 0.64 x 10~ 4 


< 2.1 x 10~ 4 


Mo = 


-250 GeV 


< 6.8 x 10" 


-7 


< 1.1 x 10~ 4 


< 8.8 x 10~ 4 


Mo = 


-500 GeV 


< 11 x 10 _ 


7 


< 0.78 x 10~ 4 


< 19 x 10" 4 


(ii) M, = 


= \M 2 = 500 GeV 


< 9.3 x 10" 


-7 


< 3.7 x 10~ 4 


< 16 x 10" 4 


(iii) ml 


= 20000 x diag{l, 1, 1, 1} GeV 2 


< 2.2 x 10" 


-7 


< 1.3 x 10~ 4 


< 18 x 10" 4 


~ 2 

m L = 


= diag{20000, 1000 2 , 1000 2 , 1000 2 } GeV 2 


< 2.4 x 10" 


-6 


< 840 x 10" 4 


< 44 x 10" 4 


(iv) ml 


= 20000 x diag{l, 1,1} GeV 2 


< 1.5 x 10" 


-7 


< 0.55 x 10~ 4 


< 2.4 x 10" 4 


ml = 


= diag{1000 2 , 1000 2 , 1000 2 } GeV 2 


< 1.5 x 10" 


-7 


< 0.59 x 10~ 4 


< 2.5 x 10" 4 


(v) ml -- 


= diag{20000, 500 2 , 500 2 , 500 2 } GeV 2 










ml -- 


= diag{500 2 ,500 2 ,500 2 } GeV 2 










Mo 


= 100 GeV, tan (3 = 2 


< 0.87 x 10~ 6 


< 500 x 10" 4 


< 51 x 10~ 4 




tan (3 = 10 


< 1.1 x 10 


-6 


< 67 x 10~ 4 


< 15 x 10" 4 




tan ,3 = 50 


< 1.2 x 10" 


-6 


< 13 x 10~ 4 


< 3.4 x 10" 4 


Mo 


= 250 GeV, tan (3 = 2 


< 2.2 x 10" 


-6 


< 450 x 10~ 4 


< 150 x 10~ 4 




tan /3 = 10 


< 3.1 x 10" 


-6 


< 88 x 10" 4 


< 50 x 10- 4 




tan (3 = 50 


< 3.4 x 10" 


-6 


< 18 x 10~ 4 


< 12 x 10" 4 


Mo 


= 500 GeV, tan [3 = 2 


< 5.1 x 10" 


-6 


< 660 x 10" 4 


< 450 x 10~ 4 




tan/3 = 10 


< 8.2 x 10" 


-6 


< 140 x 10" 4 


< 177 x 10~ 4 




tan 13 = 50 


< 9.8 x 10" 


-6 


< 29 x 10~ 4 


< 46 x 10" 4 
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FIGURES 




FIG. 1. The R-parity violating chargino-like loop diagram. 




FIG. 2. Diagram with charged gaugino/higgsino mixing. 
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FIG. 3. A diagram involving Majorana-like scalar mass insertion. 3a/ The diagram. 3b/ 
The seesaw origin of Majorana-like scalar masses for the "sneutrinos" explicitly illustrated. 
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FIG. 4. Illustrative diagrams with chirality flip on the external muon line. 4a/. A charged 
gaugino loop. 4b/ A g — A neutralino-like loop. 
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FIG. 5. Contours of B([i — > cy) in the (real) plane of (ji 3 , A 321 ). The 90% C.L. allowed region 
is shaded. 
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FIG. 6. Contours of B(fj, — ► e7) in the (real) plane of (fj, x ,iJ, 2 ). The 90% C.L. allowed region 
is shaded. Note that the approximation we used for the external lepton lines is less applicable at 
the right and top ends of the plot where the result should be read with caution. 
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FIG. 7. Contours of B(p, — > erf in the (real) plane of (B 3 , A 321 ). The 90% C.L. allowed region 
is shaded. 
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Contours of B(fj, — ► erf in the (real) plane of (Bi,fj, 2 ). The 90% C.L. allowed region 
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